Nucleation in the chiral transition with an inhomogeneous background 
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We consider an approximation procedure to evaluate the finite-temperature one-loop fermionic density in the 
presence of a chiral background field which systematically incorporates effects from inhomogeneities in the 
chiral field through a derivative expansion. Modifications in the effective potential and their consequences for 
the bubble nucleation process are discussed. 

PACS numbers: 

Keyword: 



o 
o 

(N 

Vh _ 
Ph. 

< 

in 

(N 



> 

o 

(N 

O 
O 

Oh! 



'i 



The mechanism of chiral symmetry breaking and the study 
of the dynamics of phase conversion after a temperature- 
driven chiral transition can be done conveniently within low- 
energy effective models |l,2,3,4^^^i^^]. In particular, 
to study the mechanisms of bubble nucleation and spinodal 
decomposition in a hot expanding plasma flO], it is common 
to adopt the linear a-model coupled to quarks 1 1 1 ] . The gas of 
quarks provides a thermal bath in which the long-wavelength 
modes of the chiral field evolve, and the latter plays the role 
of an order parameter in a Landau- Ginzburg approach to the 
description of the chiral phase transition. The standard pro- 
cedure is then integrating over the fermionic degrees of free- 
dom, using a classical approximation for the chiral field, to 
obtain a formal expression for the thermodynamic potential 
from which one can obtain all the physical quantities of inter- 
est. 

To compute correlation functions and thermodynamic 
quantities, one has to evaluate the fermionic determinant that 
results from the functional integration over the quark fields 
within some approximation scheme. Usually one performs a 
standard one-loop calculation assuming a homogeneous and 
static background field 1 12]. Nevertheless, for a system that is 
in the process of phase conversion after a chiral transition, one 
expects inhomogeneities in the chiral field configuration due 
to fluctuations to play a major role in driving the system to 
the true ground state. In fact, for high-energy heavy ion col- 
lisions, hydrodynamical studies have shown that significant 
density inhomogeneities may develop dynamically when the 
chiral transition to the broken symmetry phase takes place 1 6], 
and inhomogeneities generated during the late stages of the 
nonequilibrium evolution of the order parameter might leave 
imprints on the final spatial distributions and even on the inte- 
grated, inclusive abundances 1 13]. 

In this paper we briefly describe an approximation proce- 
dure to evaluate the finite-temperature fermionic density in the 
presence of a chiral background field, which systematically 
incorporates effects from inhomogeneities in the bosonic field 
through a gradient expansion. (For details, see Ref. |14].) 
The method is valid for the case in which the chiral field 
varies smoothly, and allows one to extract information from 
its long-wavelength behavior, incorporating corrections order 
by order in the derivatives of the field. For simplicity, we 
ignore corrections coming from bosonic fluctuations in what 
follows. Those would result in a bosonic determinant correc- 
tion to the effective potential ifisll . To focus on the effect of 



an inhomogeneous background field on the fermionic density, 
we treated the scalar field essentially as a "heavy" (classical) 
field, whereas fermions were assumed to be "light". 

Let us consider a scalar field (|) coupled to fermions \|/ ac- 
cording to the Lagrangian 



L 



■■^\ifd,^^f-Mm^^- 



d,^^^-V{^), 



(1) 



where ^ is the fermionic chemical potential, M{(^) is the ef- 
fective mass of the fermions and y((|)) is a self-interaction po- 
tential for the bosonic field. 

In the case of the linear a-model coupled to quarks, (|) rep- 
resents the G direction of the chiral field <I> = (C),S), where TT^ 
are pseudoscalar fields playing the role of the pions, which 
we drop here for simplicity. The field \|/ plays the role of 
the constituent-quark field q= (u^d), and in = jHq is the quark 
chemical potential. The "effective mass" is given by M{(^) = 
g|a|, and V{^) = (X?/4){G^^Tf - v^f - hqC is the self- 
interaction potential for <I>. The parameters above are chosen 
such that chiral SUl{2) <S) SUr{2) symmetry is spontaneously 
broken in the vacuum. The vacuum expectation values of the 
condensates are (a) =fn and (S) = 0, where /j^ = 93 MeV is 
the pion decay constant. The explicit symmetry breaking term 
is due to the finite current-quark masses and is determined by 
the PC AC relation, giving hq = /nm^, where mj^ = 138 MeV 
is the pion mass. This yields v^ = /^ — m^fX?. The value 
of X^ = 20 leads to a a-mass, m^ = 2X^f^ + m^, equal to 
600 MeV. 

The Euler-Lagrange equation for static chiral field config- 
urations contains a term which represents the fermionic den- 
sity: 



dV 



(2) 



and the density of fermions at a given point xq has the form 



p(xo) =Sp{xo 



G-e' 



-M{x) 



XQ 



(3) 



where \xo) is a position eigenstate with eigenvalue xq, and Sp 
represents a trace over fermionic degrees of freedom, such as 
color, spin and isospin. 

In order to take into account inhomogeneity effects of the 
chiral background field, (|), encoded in the position depen- 



dence of M in 01, we resort to a derivative expansion as ex- where one should notice thatxo is a c-number, not an operator, 
plained below |14]. 



In momentum representation, the expression for the 
fermionic density assumes the form 



^'*'=^''^?/(w«-' 



1 



f{mn^!^)-l-k^M{x) 



Jk-x 



(4) 



One can transfer the xq dependence to M{x) through a unitary 
transformation, obtaining 



d^k 



^(«)=S,Tl,l^ 



1 



{2n)^ f(^ia)^^jii)-y.k^M{x^xo) 



Now we expand M{x-\-xo) around xq : 



M(x + xo) = M(xo)+AM(x,fo) 



M{xo) + ViM f + -ViVjM x'xJ + • • • , (6) 



(5) and use x^ = —i^kt to write 
J 



d^k 



"^'"^'"^lini 



{27zy f{ia)^^lu)-y-k^M{xo) 



1 -1 



l^AM{-Nk.,xo) 



f{mn^l^)-y-k^M{xQ) 



(7) 



r 



To study the dynamics of phase conversion after a chiral small ammount of momentum to the chiral field. Under this 

transition, one can focus on the long -wavelength properties assumption, we can expand the expression inside brackets in 

of the chiral field. From now on we assume that the static Eq. in a power series: 
background, M(x), varies smoothly and fermions transfer a 



pW-^o^E/ll 



(2%f f{m^ +^) -y-k^M{x) 



z: U-^Y 



AM(-/V;^.,f) 



f{mn^ln)-y-k^M{x) 



(8) 



Eq. (|8|), together with 

AM(-/V^^.,f)=V,-M0') V^^^. + iv.V.-MQ') Vkykj^'-- , 

(9) 
provides a systematic procedure to incorporate corrections 
brought about by inhomogeneities in the chiral field to the 
quark density, so that one can calculate p {x) = po (•^) + p i (•^) + 

p2(-^) H order by order in powers of the derivative of the 

background, M{x). 

The leading-order term in this gradient expansion for p{x) 
can be calculated in the standard fashion fl21 and yields the 
well-known mean field result for the scalar quark density 



Po 



./ 



d^k M{(^)/Ek{(\>) 



- {f^q -^ -l^q) 



(10) 



where v^ = 12 is the color- spin-isospin degeneracy factor, 
Ek{^) = (P +m2(^))1/^ and M(4>) = g\(^\ plays the role of 
an effective mass for the quarks. The net effect of this leading 
term is correcting the potential for the chiral field, so that we 



can rewrite Eq. (O as 



V^C^: 



av. 



eff 



d^ 



(11) 



where Ve/f = y(^) + Vq{(^) and 



V'^^-^'^Tl-0^ln{e^^^^^y^^y^ + l)- 



(12) 

This kind of effective potential is commonly used as the 

coarse-grained thermodynamic potential in a phenomenologi- 

cal description of the chiral transition for an expanding quark- 

gluon plasma (3111113]. 

The next non-trivial term in the expansion contains two 
contributions: one coming from pi and another from p2. This 
is due to the rearrangement of powers of the gradient oper- 
ator. This term will correct the Laplacian piece in the chi- 
ral field equation. Dropping zero-temperature contributions 
which can be absorbed by a redefinition of the bare param- 
eters in y, a long but straightforward calculation yields 114,1 
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FIG. 1: Yeff and Yin% for different values of the temperature T = [108 (a), 116 (b), 124 (c), 132 (d)] MeV at ^/^ = and for g = 5.5. 



{pi+P2) = -{V^M)W,{T,n„<^) 



(13) 



where 



W,iT,^„<\,) = ^f dkk^[H{Ek,T,iJ,)- 



-H{Ek,T,-iJg)], 

(14) 

and H{Eii,T,iJq) is a complicated function whose particular 
form is not very illuminating. In the low-temperature limit 
(PM >> 1), the integral above is strongly suppressed for high 
values of k, and the leading term has the much simpler form 



Wg{T,p„<^) 



2% 



-PI^I(p|M|)3/2, 



(15) 



which gives a a better idea of the profile of the first inhomo- 
geneity correction. One can already anticipate that it will be 
concentrated in the same region where the homogeneous cor- 
rection was significant, i.e. PM < 1 (cf. ||5J), being expo- 
nentially suppressed for higher values of the field. In fact, a 
numerical study of the complete Wq shows that this function 
is peaked around (|) = and non-negligible for P|(|)| < 1 1 14]. 



The Euler-Lagrange equation for the chiral field up to this 
order in the gradient expansion reads 



V^4>: 



[l^gWq{T,lAqM 



'""'^'^^^^ ^vlM • m 



The complete new effective potential can be obtained from 
our previous results by numerical integration. In order to pro- 
ceed analytically, though, we choose to fit its derivative, which 
we know exactly up to this order, by a polynomial of the fifth 
degree. Actually, we know that the commonly used effective 
potential, Veff = V + V^, can hardly be distinguished from a fit 
with a polynomial of sixth degree in the region of interest for 
nucleation |5, 16]. Working with fits will be most convenient 
for using well-known results in the thin-wall approximation 
to estimate physical quantities that are relevant for nucleation, 
such as the surface tension and the free energy of the critical 
bubble 1 17]. 

We can now integrate analytically the polynomial approxi- 
mation to the derivative of the complete effective potential. In 
Fig. 1 we display the curves for Veff and Vtn for a few values 
of temperature and jUq = 0. From the figure one can notice 
a few consequences of the inhomogeneity correction. A first 



general effect is the smoothening of the effective potential. In 
particular, and most importantly, the barrier between the sym- 
metric phase and the broken phase is significantly diminished, 
as well as the depth of the broken phase minimum. Therefore, 
one can expect an augmentation in the bubble nucleation rate. 
In principle, one should also have better results from calcu- 
lations within the thin- wall approximation. Also, the critical 
temperature moves up slightly. 

Let us now consider the effects of the first inhomogene- 
ity correction on the process of phase conversion driven by 
the nucleation of bubbles BITIi . To work with approximate 
analytic formulas, we follow Ref. |5] and express Vtn over 
the range < (|) < T in the familiar Landau- Ginzburg form 
y&ff ^ Ln=o ^n ^" . Although this approximation is obviously 
incapable of reproducing all three minima of Vin, this polyno- 
mial form is found to provide a good quantitative description 
of Vin in the region of interest for nucleation, i.e. where the 
minima for the symmetric and broken phases, as well as the 
barrier between them, are located. Using the thin- wall approx- 
imation, we can compute the relevant quantities in the way 
described in Ref. |5]. Nevertheless, it was shown in |5] that 
the thin- wall limit becomes very imprecise as one approaches 
the spinodal. In this vein, the analysis presented below is to be 
regarded as semi-quantitative. To be consistent will compare 
results from the homogeneous calculation to those including 
the inhomogeneity correction within the same approximation 

Q. 

To illustrate the effect, we compute the critical radius, R^ 
and free energy of the critical bubble in units of tempera- 
ture, Fh/T, for three different values of the temperature. For 
r = 108 MeV, corresponding to the spinodal, which is not 
modified by the first inhomogeneity correction, the corrected 



values are Re ~ 0.98 fm and Fh/T ^ 0.20, as compared to 
Re ^ 1. 1 fm and F^/T ^ 0.9 in the homogeneous case. The 
same computation for T = 116 MeV yields Re ^2.15 fm and 
Fb/T ?^ 1.14, as compared to Re ^ 2.2 fm and Ft/T ?^ 2.1. At 
r = 123 MeV, which corresponds to the critical temperature 
for the homogeneous case, the critical radius and F^/T diverge 
in the homogeneous computation, whereas 7?^ ^ 35 fm and 
Fij/T ^ 394 including inhomogeneities. The numbers above 
clearly indicate that the formation of critical bubbles is much 
less suppressed in the scenario with inhomogeneities, which 
will in principle accelerate the phase conversion process after 
the chiral transition. 

In summary, the first inhomogeneity correction brings sig- 
nificant modifications to the form of the effective potential. 
Besides a general smoothening of the potential, the critical 
temperature moves upward and the hight of the barrier sepa- 
rating the symmetric and the broken phase vacua diminishes 
appreciably. As a direct consequence, the radius of the criti- 
cal bubble goes down, as well as its free energy, and the pro- 
cess of nucleation is facilitated. This effect will compete with 
dissipation corrections that tend to slow down the process of 
phase conversion 1 16] . Although the numbers presented above 
should be regarded as simple estimates, since they rely on a 
number of approximations, the qualitative behavior is clear. 
A more detailed analysis can be found in Ref. 114] 
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